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Abstract This report is a review of Darwin's classical theory of bodily tides in which we present 
the analytical expressions for the orbital and rotational evolution of the bodies and for the energy 
dissipation rates due to their tidal interaction. General formulas are given which do not depend on any 
assumption linking the tidal lags to the frequencies of the corresponding tidal waves (except that equal 
frequency harmonics are assumed to span equal lags). Emphasis is given to the cases of companions 
having reached one of the two possible final states: (1) the super-synchronous stationary rotation 
resulting from the vanishing of the average tidal torque; (2) the capture into a 1:1 spin-orbit resonance 
(true synchronization). In these cases, the energy dissipation is controlled by the tidal harmonic with 
period equal to the orbital period (instead of the semi-diurnal tide) and the singularity due to the 
vanishing of the geometric phase lag does not exist. It is also shown that the true synchronization with 
non-zero eccentricity is only possible if an extra torque exists opposite to the tidal torque. The theory 
is developed assuming that this additional torque is produced by an equatorial permanent asymmetry 
in the companion. The results are model-dependent and the theory is developed only to the second 
degree in eccentricity and inclination (obliquity). It can easily be extended to higher orders, but formal 
accuracy will not be a real improvement as long as the physics of the processes leading to tidal lags is 
not better known. 
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synchronous rotation, stationary rotation, capture into 1:1 resonance 
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1 Introduction 

This report aims at presenting the main ideas of Darwin's classical theory of bodily tides (Darwin, 1879, 
1880) and its consequences in the study of tidal friction effects on close-in satellites and exoplanets 
(hot Jupiters). In that sense, it is not original: the works of Goldreich (1963), Kaula (1964), Alexander 
(1973), Zahn (1977), Mignard (1979), Hut (1981), Eggleton et al. (1998), etc. already explored the 
consequences of Darwin's theory. They showed that Darwin's theory is sufficient to understand the 
main effects of tidal friction in the Solar System. Indeed, this report was initially written as an attempt 
to have a document presenting the fundamental equations of tidal friction in a simpler way, close to the 
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approach followed by Jeffreys (1961). Kaula's very complete theory with consideration of higher order 
tides, and with many infinite series in eccentricity and inclination, must be handled with enormous 
care to avoid being lost in the successive summations and in the interpretation of some coefficients not 
unambiguously defined. An additional difficulty arises from the fact that many formulas appearing in 
the literature citing Kaula (1964) are not actually found in that paper. They may have been derived 
from those in the paper, but implicitly using additional assumptions (on lags, dissipation and rotation 
frequencies, for instance). These additional assumptions are never mentioned in these indirect citations. 
More simple and self-contained approaches to the problem are certainly useful. They are found in the 
literature but, most of them consider only parts of the subject and the existing results are fragmented 
in a large number of different papers. 

The main differences of this report with respect to Darwin's work concern the phase lags of the 
tidal waves. In Darwin's theory, each phase lag is assumed to be proportional to the corresponding 
wave frequency. Here, the main equations are obtained considering each lag as an independent quantity, 
only assuming that equal frequencies lead to equal lags and that the lag vanishes when the frequency 
goes to zero. Results with assumptions on the lags are introduced afterwards, almost at the end of the 
paper, thus avoiding the need for reworking the whole theory when a different assumption is adopted. 
The equations given here for the orbit evolution do not depend on particular assumptions on the 
lags' dependence on frequencies. Another important difference comes from the fact that Darwin (1880) 
considered simultaneously with the tidal effects those arising from the oblateness of the deformed body. 
In his theory, the orbital plane precesses due to the oblateness. 

A major difference between this report and the traditional literature is the consideration, here, of the 
two different final states for the rotation of close-in companions. The first, called "stationary rotation", 
corresponds to the equilibrium situation reached by the close-in companion in which the average tidal 
torque acting on it vanishes and there is no other torque acting on the body. As it is well-known, 
if the orbital eccentricity of the companion is not zero, the resulting stationary state is a rotation 
slightly faster than the orbital motion (super-synchronous rotation). This is the situation adopted in 
the majority of studies of tidal evolution of planetary satellites and is expected to happen in the case of 
a fluid companion. The second, called "synchronous rotation" , corresponds to the equilibrium situation 
reached when the companion has an important permanent, solid-body-like, equatorial bi-axiality ( J22 
or Jai-component of its gravitational potential). In this case, tidal friction drives the rotation close to 
synchronous rotation allowing for the companion to be captured into a 1:1 spin-orbit resonanceQ due 
to the torques arising from the equatorial bi-axiality. This is the situation adopted in the majority of 
studies of satellite rotation (Cassinian theories). It can be achieved by planetary satellites and close-in 
super-Earths. 

In the study of orbital evolution, we preferred an approach different from that generally followed. 
Instead of using the direct Gauss or Lagrange equations to obtain the variation of the orbital elements, 
we preferred to obtain these variations from the application of conservation laws. The derivation is not 
as direct as with Gauss or Lagrange equations, but we consider it as important to keep in touch with 
the underlying physics. As a bonus, the equations giving the dissipation of energy in the interior of 
the deformed bodies are obtained as part of the process. 

In addition, we mention that some studies consider only the force acting on one body due to the 
tidal deformation of the other and neglect the bulk reaction force acting on the deformed body (they 
consider only the torque due to that reaction). In most of such theories, this is not an error: They 
aim at studying the friction due to tides raised by a satellite on a planet, and the acceleration due to 
the reaction force has to be divided by the mass ratio planet-to-satellite. The effects are thus much 
smaller than the direct ones, but the results obtained with this assumption cannot be easily applied 
to different cases. 

Finally, it is also worthwhile mentioning that, following Darwin, we do not use Love theory (which 
did not exist when his works were published); it is noteworthy that the results usually obtained with 
Love theory do not essentially differ from those obtained just by adopting the first-order Jeans spheroid 
as equilibrium figure of the body under the action of an external attraction combined with Darwin's 
ideas on the lag and reduced height of tidal waves due to the visco-elastic nature of the bodies. One 
important difference of the approach followed in this report with respect to many other works on 
Darwin's theory is the possibility of introducing different coefficients for harmonics of the tidal wave 

^ We should add the possibility of capture into a different spin-orbit resonance long before approaching 
synchronization, as in the case of Mercury (see Lemaitre et al. 2006). 
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Fig. 1 Equilibrium figure of tlie deformed body m under the gravitational attraction of M 

with different frequencies, instead of introducing the same "Love number" for all harmonics arising 
from U2 (see eqn. I14[) . Anyway, the coefficients used in the theory are of the same nature as Love 
numbers and Love theory (Love, 1927) should be used to estimate their values (but such estimation is 
not done here). 

This report is divided into many small sections. Sections 2 to 7 include the necessary derivation 
of the forces between the two bodies due to the tide raised on one of them by the action of the 
other. In two of these sections (sections 3 and 4), the different tidal waves are described and the rules 
used to introduce tidal wave lags are presented. "Synchronous" and "stationary" rotation of close-in 
companions are discussed in Sections 8 and 9. Sections 10 to 13 present the effects due to tides raised 
on close-in companions. Section 14 presents the effects of tides raised on the central body. In section 
15, the results are used to reproduce results from linear theories with a constant time lag. At the end, 
sections 16 to 18 discuss how the two tides must be added to obtain the variation of the orbital elements 
of the system. For sake of comparison to other works, section 16 gives the variation of the mean-motion 
and eccentricity under different hypotheses, with emphasis on linear theories with a constant time lag. 

2 The static equilibrium tide 

We consider two bodies orbiting one around another, separated by a distance r, and assume that one of 
the bodies is deformed under the tidal action of the other. In the whole report, m is the deformed body 
and M is the outer mass (a mass point) responsible for its deformatior0. We choose a reference frame 
with origin in the tidally deformed body and set r and r* as the position vectors of, respectively, M 
and an arbitrary point M* in space. The body m is initially considered as a homogeneous and perfect 
inviscid fluid that assumes the equilibrium shape dictated by its internal gravity and the tide generating 
potential due to the external mass Al. If rotation is neglected, the equilibrium figure is a Jeans spheroid 
(c — b) whose axis of symmetry is pointed towards M (see figure [1} and whose prolateness is 



(Tisserand, 1891). R is the mean radius of m. If the body rotates and the external mass lies on the 
body's equatorial plane, the equilibrium figure changes into a Roche ellipsoid the smaller axis of which 
is directed along the rotation axis (see Chandrasekhar, 1969). The prolateness of the equator is the 
same as given above and it is, generally, almost invariant to the rotation velocity (see Appendix A). 
The forces arising from the polar oblateness may be, in a first approximation, superimposed on the 
tidal forces. 

The potential raised by a prolate spheroid of mass m at an external point located in r* is given by 



^ When establishing the main equations, we do not need to identify which is which in the system of two 
bodies. The relative size of the bodies may be considered only in the applications because the equations are 
valid even when the bodies have comparable masses. To take into account the tides raised on both bodies, it 
is necessary to add the two effects following the equations given in the three last sections of the report. 




(1) 




(2) 
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(see Beutler, 2005) where A, B, C are the moments of inertia of m with respect to its principal axes. 
{A is the moment of inertia with respect to the symmetry axis of the spheroid; C = B.) We note that, 
since the symmetry axis of the spheroid coincides with the direction of r, •Z' is the angle formed by the 
position vectors r and r*. Assuming that B is proportional to + (? and A is proportional to &^ + c^, 
there follows B — A :^ eA. Hence, 

— (3cos^«^-l) (3) 
or, introducing the parametei0 kf — IbA/AmR^ , 

u = -- ^^3;:;:^(3cos2«'-i); (4) 

In order to express cos If" as a function of the components of the vectors r, r*, we choose a spherical 
coordinate system (figure [2]) so that r = {r,6,ip) and r* = {r* ,9* ,ip*). The two angles considered for 
each point are their co-latitudes {0,9*} and longitudes {(p,(p*). We have 

cos 'F — cos 9 cos 9* + sin 9 sin 9* cos{tf — f*)- (5) 

Let us now introduce the relations between r, 9, (p and the orbital elements of M. From the spherical 
triangle shown in the lower part of figure [21 we have 

sin 9 = cos(cij + v) cos ip + sin(cij + v) sin tp cos / (6) 



cos 9 — sin(w + v) sin / (7) 

and 

(^~t) + Lj- -sin(2t; + 2cj)sin^/ + 0(/'') (8) 

where lu is the argument of the periapsis of M and v is the true anomaly. / is the obliquity, that is, 
the inclination of the orbital plane with respect to the reference plane (equator) . From the equations 

^ This parameter is often called fluid Love number because it is equal to the secular Love number of a rotating 
Maclaurin spheroid whose oblateness is e. For a homogeneous sphere, kf — 1.5 (since, in this case, the moment 
of inertia is A = OAmR^.) 
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of the Keplerian motion we have, in terms of the mean anomaly i, to second order in the eccentricity 
e: 



V = ^ + 2esin£ + -e^ sin2€, 
4 



1 - ecos^ + ie^(l - cos2^) 



(9) 



(10) 



The next step is to substitute the above equations into ^ and ^ and expand to O(e^) and to 0{P). 
For sake of simplicity we adopt the notations 



S = sin/ 
P = sinr 
g = sin 20* 

At the orders considered, the static equilibrium tidal potential at the point r* is 



(11) 
(12) 

(13) 



t/2 = 



■ikfGMW 



1 - ^e^ - ^52^ P2 cos(2(y5* -2t- 2uj) + ^eP^ cos(2(^* - U - 2lo) 

1 17 

--eP2 cos(2</j* - £ - 2cj) + ye^P^ cos(2(p* - 4£ - 2uj) 

2 - 3P2^ e cos £ - (^3 - ^P2^ cos 2£ 

+gS'( sin V?* -sin((^* -2£-2w)) + ip^S^ ( cos 2^?* + cos(2^ + 2a;) 



(14) 



This completes the computing of the potential raised by the tidally deformed body on an arbitrary 
point in space up to the second order in eccentricities and inclinations. 



3 Tidal waves 

To interpret the terms in t/2, we consider a point fixed to the surface of the body. The longitude of this 
point is Lp* — Qt + ifQ, where CI = 12 k is the rotation angular velocity vector of m assumed normal 
to the reference plane and <Pq is a constant. Each term depending on ip* in C/2 corresponds to a tidal 
wave traveling on the body with given direction and velocity. All terms contribute, in different ways, 
to the formation and evolution of the tidal bulge on the body. The terms of U2 may be divided into 
three groups: 

1. - Sectorial components having the form P^ cos{2ip* — a), where a is the corresponding phase. 
The amplitude of these terms is maximum when P^ is maximum, i.e. at the equator, and at the 
longitudes ip* — a/2 and is decreasing towards the poles. The main term (No. in Table [T]) is a 
wave with period 7r/(i7 — n) (i.e. half the synodic rotation period) with two maxima located on the 
intersection with the equator of the meridian passing through the sub-M point and the other on its 
antipodal. If n <C the period is nearly half of the rotation period. The next term (No. 1 in Table 
[T]) is a wave with period 7r/(i7 — 1.5n) (i.e., larger than half synodic rotation period). The wave 
has two maxima which are located on two antipodal points on the equator. One of them lies on the 
meridian of the sub-M point when the tide generating body is at the periapsis i — and the other 
when i — n. The high tide moves, in this case, more slowly than the sub-M point. Similar analyses 
can be done for the other terms showing the argument 2ip*. They are often called semi-diurnal 
tides since, on the Earth, they have periods close to 12 hours. They are shown in Table [U which 
also summarizes the interpretation to be given in other cases {Q <^n and I? ~ n). 
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No. 


frequency 


Type 1 

O n 
J i 1 1 


Type 2 

J Ci ( I 


Type 3 


1 7\ 1 




1 TTTT 7% 1 


1 ' — p 1 1 

semi-diurnal 




1 ' 1 1 

semi-annual 


1 


2n - 3n 


semi-diurnal 


monthly 


3"* of annual 


2 


2n - n 


semi-diurnal 


monthly 


annual 


3 


2n - An 


semi-diurnal 


semi-monthly 


4*" of annual 


4 


2n 


semi-diurnal 


semi-monthly 


"semi-diurnal" 


5 


n 


monthly 
(radial) 


monthly 
(radial) 


annual 
(radial) 


6, 7 


2n 


semi-monthly 
(radial) 


semi-monthly 
(radial) 


semi-annual 
(radial) 


8 


n - 2n 


diurnal 


monthly 


semi-annual 


9 




diurnal 


monthly 


"diurnal" 



Table 1 Tidal potential analyzed term by term. The given tide frequencies and corresponding names refer 
to how the tidal potential is felt on a given (fixed) point of the body. For the type 2 tides, the paradigm is 
the Moon, but when the synchronous companion is an exoplanet, the names annual and semi-annual would be 
more appropriate. 



2. - Zonal components independent of the longitude ip* . These components of the tidal potential 
oscillate all over the body with amplitudes depending on the latitude of the points and on the mean 
longitude of the tide generating body. For example, the component proportional to (— 2-|-3P^) cos£ 
has relative maximum amplitudes at the equator and the poles, but with inverted phases. The phase 
inversion occurs at the critical latitude corresponding to — 2/3 {9* — 54.7 degrees). These terms 
are often called radial tides because there is no propagation of a crest around the body. 

3. - Tesseral components having the form QSsv[i{lp* — a). The amplitude of these terms is maximum 
when Q = 1, i.e. at the latitudes ±45 degrees {9* — 45 and 135 degrees) and at the longitudes 
(/9* = a, with a phase inversion on the equator. They are often called diurnal tides since, on the 
Earth, they have periods close to 24 hours. 



In table [U we summarized the analysis of the tidal waves for three different cases depending on 
the rotation speed of the deformed body. Type 1 corresponds to a body rotating with angular velocity 
much larger than the orbital mean motion (i? 3> n). It is the case of the Earth-Moon system, with the 
Earth as the deformed body and the Moon as the perturbing one. Type 2 corresponds to synchronous 
or almost synchronous motions and, again, the Earth-Moon system serves as an example, but now 
the Moon is the deformed body and the Earth is generating the tide. Looking at Table [H we see that 
synchronization gives rise to terms whose period is related to the rotation period of the companion, 
and they are called monthly, semi-monthly, etc. The names come from the tidal action of the Earth on 
the Moon (for this reason, the semi-monthly tide is often called fortnightly). Since fl n,we could also 
use the names diurnal, semi-diurnal, etc. (referring to the companion's rotation). Type 3 corresponds 
to a slow rotating body {Q <^ n)] it is the case of the tides on a typical main sequence star due to 
a close-in planet (hot Jupiter). Using names similar to those used in the other cases and taking into 
account that the main period is the planet's orbital period (or "year"), we will call them, respectively, 
annual, semi-annual, tierce-annual and so orQ. 

It is worth emphasizing that the given tidal frequencies and corresponding names refer to how 
the tidal potential is felt on a given (fixed) point of the body. The propagation of the tidal wave in 
the body must be analyzed separately. For instance, on the Earth, the tidal bulges of both diurnal 
and semi-diurnal tides circulate around the Earth with the synodic rotation speed. The names and 
frequencies given in Table [T] refer rather to the shape of the tidal wave. 

Table [T] is limited to the tidal components appearing in the given expansion of U2- When higher 
order terms are considered, many other frequencies appear. 



* We have, however, to keep in mind that these "years" are very short. For instance, the orbital period of 
OGLE-TR-56b, one of the shortest known, is only 1.21 days. In type 3 tides, "diurnal" is much slower than 
" annual" . 
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4 The tidal phase lags 

In the previous sections, we have considered that m is a perfect inviscid fluid that reaches the equihb- 
rium figure instantaneously under the attraction of M. However, in a real body, the viscosity introduces 
a delay between the tidal action and the corresponding response. 

The main characteristic of Darwin's theory is to consider the potential U2 as a composition of 
periodic terms with different frequencies and to introduce in each periodic term a delay in the form of 
a lag angle (Darwin, 1880). U2 is then expanded to first order in the lags. The trigonometric functions 
are expanded in the following way: 

cos(^i - £i) ~ cos^i + El sin^j (15) 
sm{$i — Ei) ~ sin<?i — eiC0S<?i (16) 

where is a generic time-dependent argument. 

This is not the only way of introducing the lags. We may mention the theories of MacDonald 
(1964) and Mignard (1979, 1980) as paradigms of different approaches. In both cases, the tidal lag is 
associated with the displacement of the tidal bulge dragged by the rotation of the body. In MacDonald's 
theory, the tidal lag is a constant and, in Mignard's theory, the tidal lag is proportional to the relative 
(synodic) rotation speed. Eggleton et al. (1998) use a different physical approach which, however, 
leads to the same tidal force given by Mignard (1979). When these theories are used to study Earth's 
bodily tides, their results do not differ essentially from those obtained by Darwin. Indeed, on Earth, 
the tidal effects are dominated by semi-diurnal tides, and the different approaches address more or 
less the same problem. Difficulties arise in the case of synchronous and quasi-synchronous rotation. 
The classical theories show that, in absence of non-tidal torques, no synchronization is possible if 
the orbital eccentricity is not damped to zero. The tidal torque on the deformed body vanishes for 
a rotation velocity slightly larger than the orbital mean motion. A simple physical reasoning shows 
that this result is expected. Indeed, the tidal torque is inversely proportional to the sixth power of 
the radius vector (see eqn. (55]). When the orbit is eccentric, the torque will be much greater at the 
periapsis than at the apoapsis, and, thus, the average will correspond to a rotational angular velocity 
exceeding the orbital mean motion. Therefore, the rotation velocity in the stationary solution is larger 
than the orbital mean motion. Linear theories with a constant time lag give for the stationary rotation 
velocity: 

n = n{l + 6e^) (17) 

(see eqn. l36l) . In the case of Titan, this result corresponds to i? = n(l -I- 5 x 10~^), that is, to a synodic 
rotation period of about 8.5 yrs, which would be observable notwithstanding the difficulties in the 
identification of features in the surface of Titan. Several attempts were done to modify Mignard's and 
similar theories by introducing phase lags depending nonlinearly on the wave frequencies (see e.g. Sears 
et al. 1993). 

The question of the law that should be used to relate tidal lags and frequencies is a controversial one. 
Some results obtained from the study of the Moon (Williams et al. 2006) indicate that the variation, 
for a large range of frequencies, is small: the quality factoid Q increases from 30 for one month to 34 
for one year. This corresponds to a power law with an exponent ^ 0.04. On the other hand, Efroimsky 
and Lainey (2007) collected geophysical data which indicate, in the frequency range of our concern, an 
inverse law: the dependence of lag on the frequency follows a power law with a negative exponent in 
the range (-0.4,-0.2). 

Darwin's theory introduces naturally different lags for different terms without the need of intro- 
ducing a priori a particular law (notwithstanding the fact that Darwin himself used a linear law). We 
will just assume that equal frequencies correspond to equal lags and that the lag vanishes when the 
frequency tends to zero. We also pay attention to the fact that some frequencies in Table [T] become 
negative when /? ~ n or i? <C n (Types 2 and 3 tides). The rule is that the actual phase of the terms in 
the dynamic equilibrium tide lags behind the corresponding phase in the "static" case. Thus, in type 
2 tidal interaction, the lags may be such that Si '~ —£2 < and eg '~ —eg < 0. In the same way, in 
type 3, we have < for all subscripts corresponding to negative frequencies (i = 0, 1, 2, 3, 8). 

^ The quality factor is usually introduced through its relationship to the geometric lag angle A: Q ^ 1/2 A 
(see MacDonald, 1964; Efroimsky and Lainey, 2007); the corresponding relationship when the phase lag of the 
semi-diurnal tide is used instead of the geometric phase angle is Q ~ l/eo- This definition needs to be amended 
when £0 becomes small. See sect. 111.3] 
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5 The dynamic equilibrium tide 



In addition to the phase lag of the tidal waves, we may assume that the body does not reach total 
deformation and thus substitute, in the coefficients, the factor fcy by dynamical counterparts ki which 
may be assumed to depend on the frequencies of the corresponding tidal waves. However, in first- 
order theories, every ki appears in the tidal potential multiplied by the corresponding lag Si. We will 
put in factor one value kd (called dynamic Love number) for the main tidal wave (which, for the 
Type 1 of tidal interaction may be taken as the semi-diurnal Love number k) and merge the others 
with the corresponding phase lag. The lags Si and the response factors ki are very distinct physical 
quantities but, as far as only the potential due to m is considered, they do not need to be considered 
separately. We may keep in mind that when physical interpretations are sought, it is easy to separate 
their contributions. To avoid any misinterpretation, we introduce the modified lags e'^ = {kj/kd)sj and 
use them instead of Ej. 

Thus, instead of t72, we have 

where U2 is given by eqn. (fT4| and 
SkMMR^ 



U2 = U° + f/lag 



(18) 



C/lair — 



P^e'o {2-5e-'-S^) sm{2ip* -21- 2w) 



+eP- 



' (7e[ sm{2ip* - 3i - 2uj) - e'2 sin(2(^* - £ - 2uj)^ 

+17e^P^e'^ sin(2v5* -U- 2w) + P'^S^e'^ sin 2ip* 
-eeg (4 - ) sin £ - Se^e^j (2 - SP^ ) sin 2£ 

+P^S^e^ sm{2£ + 2uj) + 2Qs{e'^ cos{ip* -2i- 2uj) - 4 



g COS f 



(19) 



We recall that, in classical theories, in which tides on the Earth are the only ones considered, the lags 
are introduced as 2e because the main terms correspond to semi-diurnal tides (see Jeffreys, 1961). In 
this report, all lags are introduced as Si, in a unified way. 



6 Tidal forces acting on the tide generating body 

The perturbing force acting on a point of mass M* placed in r* = {r* ,9* ,(p*), due to the disturbing 
potential U2 is given by 

F = -M*grad . U2 = -M*— r'* - — ^ 9* - ^ if'* (20) 

r ^ dr* ^ r* ^ 89* ^ r* sin 9* dip* 

Fi F2 F3 (21) 

where {r*,9*^ip*) form, at r*, a right-handed orthogonal set of unit vectors in the positive direction 
of the increments. (95* = r* x 9*.) The total force is decomposed into three orthogonal components. 
Knowing that gradj.. t/2 = gradj.. + grad^.C/zag, we may calculate the forces separately for each part 
of the disturbing potential and add them later. These forces are given in Appendix C of the astro-ph 
version of the report. They are important intermediaries in the calculations, but have no importance 
per se. 

To obtain the forces on the tide generating body M, it is enough to make the identification of 
(M*, r* ,9*,(p*) and (M, r, 9, (p). For sake of completeness, we note that this identification is done after 
the calculation of the gradient of U2- 

The forces derived from U2 are 

ri'.-?Mpl, F|' = F," = o. ,22, 
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Only the radial force survives in the absence of lags, as expected from the symmetry of the resulting 
configuration (see figure [1]) . 

Making the identification in the forces arising from J/iag and introducing the expressions given by 
eqns. ([5l)- (fT0)) . we obtain the force components 



Fi(r) 



8a7 



3e(84 - 7e; 



24) sin ^ 



(23) 



-3e^(21eo - Ae'^ - 17 e'^ + 44 + Ssg) sin2€ + iS'^ie'^ - 4 - 4 - 24 + 2e'Q) sm{2l + 2uj) 



F2{v) = 



8a7 



eS'(+84 - 74 - £2 + 64 - 12£8 + 44 + 645 + Se'ie) coscj 



(24) 



-eS'(84 - 7e'i - 4 + 665 + Ae'^ - 12e(, - 6e'i4 + f 4£'i7) cos(2£ + uj) - 4S'(4 - e^,) cos(£ + w) 



^3(r) 



3kdGM^R^ 
8^ 



(4 - Ue^ - SS^Wo + 56e^e[ + 2S^{e'g + 4) + e(f 64 + 14e; - 2e'2) cos£ 



+e2(44e;, - 84 + 34e^) cos2£ + S^ie^ + 24 - 24 - 2eg) cos(2£ + 2uj) 



(25) 



The terms with the lags e'l^^, e'l^, e'lQ, e'u introduced in F2 (cf. Errata in press) come from SsUiag (see 
Appendix D). When we assume the lags proportional to the frequencies, the forces given by eqns. (p3l 
-[55]) are equivalent to the second-degree expansion of the force given by Mignard (1979). 

The average values of the Fi can be easily found from the above equations. However, in the used 
spherical coordinates, the unit vectors are continuously changing and averages are of limited interest. 
We nevertheless note that the average radial component of the tidal force acting on M is zero and that, 
when / = 0, the component F2 vanishes. 



7 The tidal torque 

Because of the delay in the response to tide raising forces, the tidal bulge is not aligned with r. This 
causes the raising of a tidal torque 

= r X F (26) 

or, since r = (r, 0, 0), 

M = (0, M2, Ah) = -rF^e* + rF2^*. (27) 

We do not write explicitly the expressions for the two non-zero components of A^, as they are 
simply products of the force components by r. On the other hand, in the applications, we need the 
average torque, which needs to be calculated in a fixed reference frame. We thus transform Ai to obtain 
a new decomposition. This is easily done since Mi = and M3 is orthogonal to the meridian plane of 
M (see fig. [3]). If, for simplicity, we adopt here a system whose z-axis is oriented towards the ascending 
node N, we obtain 

Mx = M2 cos 6 cos (p — M3 sin (p 

My = M2C0sdsm(p + Mscosip (28) 
M^ = -M2sm6. 

The averages of M^^My, M^ over one orbital period are 
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Fig. 3 Projection of the torque components on the meridian plane of M (left) and on the reference plane 
(right) 




ref. plane 



Fig. 4 Meridian plane normal to the nodal line. Components of the torque on the tidally deformed body, k 
and ik are unit vectors. f2 is the angular velocity vector 



< M^> = 

<My> = ^^-g 5(4 + - 4) (29) 



< > 



SkdGM^R^ 



(ae'o + e\-20e'o + 49e'i + e^) + 2S'^{-2e'o + 4 + Eg)) . 



The above results deserve some comments: (1) The average torque is perpendicular to the line of 
nodes (perpendicular to the orbital plane if J = 0); (2) The results in a system whose axis is directed 
to a fixed point O (see fig. [1]) are easily obtained rotating the results by an angle ON (longitude of the 
ascending node) around the z-axis; 



8 The rotational angular momentum 

The angular momentum conservation is very stringent and cannot be used without the simultaneous 
analysis of all involved forces. For instance, the possible non-spherical shape of the bodies and the 
tidal friction in both bodies (see Section will move the planes considered in Fig. [51 In this section, 
we consider only the interaction between the orbit and the tidally deformed body. If we neglect all 
additional factors mentioned above, the change of angular momentum in the orbit may be compensated 
by a change in the angular momentum of the tidally deformed body of opposite direction. We may 
have £rot + £ = 0; but C = M. Hence £iot = —M. A conceptual difficulty in this calculation is the 
fact that we are using the equator as reference plane and it is no longer inertial if moving. To overcome 
this difficulty the simplest way is to consider the equator at a given fixed time inside the considered 
interval as reference plane. Alternatively, we could add the centrifugal and Coriolis forces acting on 
the bodies due to the motion of the reference plane. We adopt here the first of these two approaches. 
The second approach is discussed in Appendix E of the astro-ph version of this report. 

The rotational angular momentum of m is given by £i.ot ~ CJ7k where C is the moment of inertia 
of m with respect to the principal axis c (see fig. [ij , is the angular velocity of rotation of m and k 
is a unit vector along the principal axis c. 
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Since the average torque M is normal to the hne of nodes, we may decompose the equation Ciot — 
—A4 into two parts and study each part separately. In the meridian plane normal to the nodal line, 
the time derivative of the angular velocity vector is ^f2k = i7k + i7Jik where J is the variation of 

the inclination due to the variation of the unit vector k (a rotation). We introduced a new letter (J) to 
distinguish between this quantity and the variation of the inclination due to the motion of the orbital 
plane (/). From the given equation, we obtain Cf2 = ~Mz and CfiJ — —My. 
Then, taking the averages, 



8Ca6 



<f2> = - 



4e'o + e^i-20e'„ + 49^ + 4) + 2S^{-2e'„ + 4 + e^,) (30) 



<'^>^ IcQa^ ^(4 + 4 - 4)- (31) 
9 Rotation of Close-in Companions 

Let us now assume that the body m is a satellite or exoplanet orbiting close to a major primary body. 
Let us also assume that it has no other significant deformation besides the tidal one due to the central 
body. Eqn. (j30p shows that, in the first approximation, < /2 >oc Eg, indicating that the system will 
smoothly evolve towards synchronization. However, the final state is not necessarily synchronous. 

Close to synchronization, we have a type 2 tide and several lags correspond to monthly tides (see 
table [1]). We assume that they are equal, just taking as negative those corresponding to negative 
frequencies: 

4 ^ -e[ ^e',^ -e', > 0. (32) 

Hence, 



<f2> = - 



2Ca6 

\k 



(e'^-e\be'^ + l2e'^)-S\'^) (33) 



<'^>= .^o.R (4-24). (34) 



9.1 Stationary rotation 



We say that the system reaches a state of stationary rotation when the average angular acceleration 
(or the average tidal torque < >) vanishes. Solving the equation < i? > = 0, we obtain, at the 
second order in e, /, 

e'o = 12eV2. (35) 
In linear theories with constant time lag (see Section I15|) , each e[ is assumed proportional to the 
frequency of the corresponding tidal wave. When these are substituted in the previous equation, we 
obtain the well-known resul10 

^tat = n(l + 6e2), (36) 
showing that, when e 7^ 0, the rotation stabilizes at a slightly super-synchronous value. Nonlinear 
theories give similar results only modifying the numerical factor multiplying (Sears et al., 1993). 

The inclination creates an effect tending to sub-synchronize the stationary rotation, but it is of 
order 0(5*^) and therefore beyond the order of approximation adopted in this report (see Levrard, 
2008). 

The second equation resulting from the angular momentum conservation gives 

• ikdSGM'^R^ ,,2, 

<J>- - ,cna^ 4(1 -6e^). (37) 

where one mixed third order term appearing due to the substitution of Eq by its value was kept to 
show one of the few instances in which terms of this kind affect the result. 

The part of the equation corresponding to the vanishing of the torque along the nodal line indicates 
that no precession of the nodes occurs when only these forces are considered. 

® When the tidal phase lag is assumed to be frequency independent (MacDonald theory), the resulting 
stationary velocity is /?stat = n(l -I- 9.5 e^) (Goldreich, 1966). 
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9.2 Spin-orbit synchronization 

The spin-orbit synchronization condition is given by f2 — n. When e'q = 0, eqn. (I33p gives 

<f2>= 2— 6^4, (38) 



which cannot vanish if e 7^ and is thus inconsistent with a synchronization. It is worth noting that 
S does not affect the above condition (at least when terms of the fourth order in e,/ are neglected). 



9.3 Synchronous asymmetric companions 

The above conclusions were derived from the assumption that m has no other significant deformation 
besides the tidal one. If the close-in companion has a permanent (solid-like) equatorial ellipticity, 
different results follow. Indeed, in such case, we have to add the contribution of the equator asymmetry 

to C/lag, e.g. 

U22 = ^ J22P22[C0S9) C0s2(ip - ip22) (39) 

where P22 is an associated Legendre function. J22 > and (^22 are the two parameters characterizing 
the asymmetry of the gravitational fielcfl. According to Goldreich (1966), synchronous rotation will 
result when J22 is larger than a critical value (see below). In such a case, the body will end up with 
spin-orbit synchronization and, as a consequence, we have, 

£0 = 0, (40) 

instead of eqn. (l35|) . 

In this case, the equations obtained in the previous section can no longer be used. The synchronous 
motion is also a stationary solution. However it is a stationary solution of one system in which the 
companion's rotation is under the simultaneous action of two torques: the tidal torque and the torque 
due to the asymmetry of the body, which allows the motion to become synchronous. 

Since the tidal friction in a synchronous companion moving in an eccentric orbit tends to accelerate 
its motion, the body will be displaced forward but, being asymmetric, this displacement will create a 
torque in the contrary direction which will compensate the tidal torque. The average torque due to 
the displacement by an angle S forward is 

<iM22).>^ 6g^^-^J/^^Bin2^ (l-^e^-l^-). (41) 

The condition < f2 >= means, now, < + (^^22)2 0, which allows the average offset angle S 
to be determined. Using the condition e'q = 0, there follows, at the order of approximation adopted in 
this paper, 

37Vf i?3 

J22 sin 25 ~ ^fcdese^ (42) 

If we know the maximum value that 6 can reach without disturbing the capture into the 1:1 spin-orbit 
resonance, we may obtain the minimum value of J22 necessary to capture into the 1:1 resonance IlThe 
numerical exploration of some examples has shown here a large influence of the periodic terms which 
may increase the critical value of J22 by several orders of magnitude. 

The torque due to U22 also contributes to the motion of the equatorial plane of the synchronous 
companion. The y-component of the momentum is, in this case, proportional to 5J22 sin 2(5 (that is, to 
Se'^e'2) and is negligible. The x-component is 

, 3GMmR^J22S 
<M22x>= 5 — ^ cos 25. 43 

^ The sign adopted for U22 is the one more currently used (see Beutler, 2005). (^22 is the longitude of one 
point on the shortest of the two equatorial axes. 

* When the tidal phase lag is assumed to be frequency independent (MacDonald theory), the resulting critical 
value of J22 is proportional to (Goldreich, 1966). However, with the averaged equations of this section, issued 
from Darwin's theory, it appears proportional to e^. 
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This term will give rise to an acceleration in the plane defined by the rotation axis and the nodal 
line. This is the classical picture in which the rotation vector describes a precession cone and is the 
counterpart, in the body, of the precession of the nodal line given by eqn. (|66p . (The total angular 
momentum of the system is invariant to the action of J722-) The approach used here is not adequate to 
study it. Moreover, the dependence of {M22)x on 5 is of higher order and, in the first approximation, 
{M22)x can be written in a form independent of 5: 

< {M22)x > ^ ■ 44 

< {M22)x > is, thus, rather due to the figure of the companion trapped into the l:f spin-orbit resonance 
and exist even when the tide is neglected (i.e. when (5 ~ 0). The resulting precession is not an indirect 
tidal perturbation. 

It is worth mentioning that if, instead of U22, we introduce 

GmR^ 

1/31 = 2 — J3iP3i(cos6')cos((^ - (^31) (45) 

where P31 is an associated Legendre function and J31 and (^31 the two parameters characterizing the 
deformation of the equator (see Beutler, 2005), the results are very similar. We do not reproduce 
details, but just say that in this case, instead of eqn. (|^^ . we obtain 

J3isin,5'^i^^fcrf4e2 (46) 
If necessary, it is easy to consider U22 and C/31 simultaneously. 



10 The Work done by the Tidal Forces 

The work done by the tidal forces in a displacement ds is given by dW — F ■ ds, or = F • v where 
V is the velocity vector. This calculation is elementary. The only necessary precaution is to take into 
account that the usual expressions for the radial and transverse components of the velocity, 

= an (^e sin i + sin 2£j 

VT = an (1 + e cose + e'^ cos 2£ - ^e'^y (47) 

are given in a reference system whose plane is the orbit plane. The radial components in this system 
and in the spherical coordinates used in this study are the same. However the transverse component 
lying on the orbital plane needs to be decomposed along the axis of the local spherical coordinates 
frame: V2 = —vt sin (3 and = vt cos (3 (see Fig. [S]). From the triangle shown in the figure, we obtain 
(by sine and cosine laws): cos/3 = sin(^/ sin{v + to) and sin/3 — S cos(p 




The work is then easily calculated by the given scalar product. The result, averaged over one orbit 



IS 



< W >tidc = 



ZnkdGAPR^ 



4e' 



147 1 

e2(-204 + —4 + -e'2 34) - AS^ie', - e',) 



(48) 
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10.1 Energy Dissipation 

The energy release may be known from the fact that the total energy variation of the system must be 
equal to zero. Besides the orbital energy, whose variation is given above, we have the rotational energy 
of the deformed body and the thermal energy dissipated in the body. The balance equation allows us 
to calculate the energy dissipated in the deformed body. 

For that sake, we have to know the work done by the tidal torque acting on the body (which is the 
opposite of the torque acting on M), 

< w >rot = < -M ■n>= cn<h> (49) 

= - gj [4^0 + e^("204 + 49e; + 4) + 2S\-2e', + 4 + 4)_ 

If the rotation of the considered body is relatively fast (i.e. ^2 > n, Eq > 0), rotational energy is trans- 
ferred from its rotation to the orbital motion while, if it is slow (i.e. f2 < n. Eg < 0), orbital energy is lost 
and partly used to accelerate its rotation. In both cases, < W >tidc+iot < and the mechanical energy 
lost by the system must be released inside the body with the average rate < ii^th >— — < >tide+iot- 

11 Energy dissipation in close-in companions 

We consider here the case of companions having reached a final rotation with I? ~ n. In the resulting 
type 2 tide, we may take into account the condition stated by eqn. ([5^ : e'^ — —£[ ~ Eg ~ — £g > 0. 
We assume also that the lag £5 of the radial tide is equal to the lag Ej of other tidal waves with the 
same period. Hence 



< W >tidc = 



e',-e^{5e'o + 19s'^)-S'{e', + s'^) (50) 



2a6 

and < W >rot— (since J7 ~ 0). 

11.1 Stationary rotation 

If the only torque acting on the body in stationary rotation is the tidal torque, we have Eq = 12e^E2 
and the above equation becomes 

<iy>tido=i {le^ + S^)e'^. (51) 

In this case, the thermal energy released inside the companion is i?th — ~Wtidc (they are equal 
when we can neglect the changes in the equilibrium rotation speed due to the variation of the mean 
motion, which is of the order of [Rjaf' .) Therefore, the dissipation is proportional to the lag £2 of the 
monthly tide (the tide harmonic whose period equals the orbital period) . 

It is worth emphasizing that this result does not depend on any hypothesis linking lags to frequencies 
(the only assumption is that equal frequencies lead to equal lags). 

11.2 Synchronous asymmetric companions 

Proceeding as above and noting that Eg — 0, we obtain, in this case, the result 

<^>„.= _?I*|M!£!(i<,.^ + S')4. (52) 

However, if e 7^ 0, a synchronous rotation cannot exist without additional forces acting in the system, 
which will also do a work. As they counteract the tidal forces, the work is expected to have a different 
sign. From p9p . we indeed obtain 

<W>22 = (l-^e -^S)^ eE2, (53) 
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where the rightmost expression was obtained introducing the value of J22 sin 18 obtained in section 

El 

When this term is added \o <W >tidc, the resuh is 

which is formally equal to that obtained for the stationary case. It may be compared to that given by 
Wisdom (2004) and Winn and Holman (2005): 

<E,,>^ ^^j^(7e^ + n (55) 

This result only differs from the given one for using the relationship h = 5kd/3 between Love number^ 
h and kd and by adopting Q — l/e2- The difference of sign comes from the fact that Eth is the energy 
released in the body. 

The same value for the energy dissipation was obtained by Segatz et al (1988) and Wisdom (2008) 
through the direct calculation of the energy released inside the companion. It is worth noting that 
the same result can be obtained without making explicit use of < M22 >■ It is sufficient to use the 
fact that the counteracting torque must be equal (and opposite) to the mean torque accelerating the 
rotation in the synchronous case and obtain the work done by multiplying the torque by the angular 
velocity of the body (see Levrard, 2008). 



11.3 On dissipation and lags 



Lags and dissipation are different aspects of the same phenomenon. Viscosity driven tidal friction is 
directly responsible for the production of heat inside the body and, at the same time, by delaying the 
response of the body to external forces. From the mathematical point of view, we may fix one of them 
and obtain the other. In this paper, we introduced the lags and computed the thermal energy released. 

The results of the previous sections show how lag and dissipation are interrelated. Equations (|48p 
and (|50|) show that, before synchronization, the loss of orbital and rotational energies due to dissipation 
in the companion are, in the first approximation, proportional to the lag Eq (through Eg). Eq is the lag 
of the tide whose period is half the synodic rotation period (semi-diurnal tide). It is half the geometric 
lag (delay of the high tide with respect to the sub-M point). In terms of Eq, the quality factor Q is 
then usually defined as Q = I/sq and is singular when Sq — > 0. This singularity needs to be explained. 
In fact, the parameter that measures dissipation is not Q but 1/Q (see Munk and MacDonald, 1960, 
chap. 3), thus the real problem is not the division by zero, but the fact that the given formula leads 
to 1/Q ^ 0. 

However, looking at eans. (P5| ) and ([SO]) . we see that the right-hand sides do not vanish when Eq — *■ 0. 
Some terms indeed vanish, but others do not and the latter will be the leading terms in the resulting 
expressions oi < W >tido and < W >rot- 

The comparison of eqns. ([54|) and ([55|l has already shown that the usually adopted quality factor 
is, in this case, Q = 1/ 62- It is no longer the inverse of the lag of the semi-diurnal tide, but the inverse 
of the lag of the tide whose period is the orbital period (monthly tide). 



12 The variation in semi-major axis and mean motion 

From E — —GmM /2a and Kepler's law, we obtain, 



3n Qn'^kdMR^ 

< n > = < a > = — 

2a Qirna^ 



44 - e2(204 - ^£'1 - \e'2 + 34) - 45^(4 - 4)] (56) 



/i = 5/0/3 in a homogeneous companion. See Munk and MacDonald (1960). 
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Fig. 6 Meridian plane normal to the nodal line. Components of the torque along the z-axis (M^), projected 
on the reference plane (My) and perpendicular to the orbit {M±). u and iu are unit vectors. The orientation 
of My in this figure corresponds to the case e'o + Eg — e'g < (see egn. I29p . Note that = 0. 



12.1 Close-in companions 

In the same way, using the expression for the work given by eqn. (|5ip , we obtain for close-in companions 
in synchronous or stationary rotation, 

If we put £9 = 1/Q and 5 = 0, we obtain the usual expression given by Peale and collaborators (see 
Sect. [HI). 

If tidal lags are assumed to be proportional to the frequencies of the corresponding tidal wave, and 
if the synchronization assumption f2 — n is introduced in Eqn. (|56p . we obtain, when S = 0, the same 
result as given by Eqn. (3) of Mardling and Lin (2004) for the variation of the semi-major axis due to 
the tides in a synchronous companion[3- However, the result thus obtained does not take into account 
the impossibility of having f2 = n without the existence of an additional non-tidal torque counteracting 
the tidal torque. The consideration of this necessary torque reduces the term 19e^ appearing in the 
equations for d and fi to 7e^. 



13 Variation in Eccentricity and Inclination 

The orbital angular momentum of the system is given by £ = I'Cj u where u is the unit vector perpen- 
dicular to the orbital plane and 

Mm o / GM 



1 „i ivi a o / ^jrivi n I — , , 

= — -T7"aVl-e2 = v/l-e2. 58 

m -I- Af na 

a, e are, respectively, the semi- major axis and eccentricity of the relative orbit (astrocentric). 

The fact that the torque and the angular momentum are both perpendicular to the line of nodes 
allows the equation £ = to be reduced to the plane defined by them. In that plane, the time 
derivative of the angular momentum is L — ■^^\C\vi -\- |yC|/m and, from the given equation, we obtain 

= A/_i and \C\I = M||, where M_l and M|| are the components of along the directions of vl 
and iu. 

From figure [SI we have 

Afx = cos / - Mj^ sin / 

Af|| = -M^_ sin / - My cos /. (59) 

Solving the two components of the equation £ = A/( , we obtain 

SnekdMR^ , , 49 , 1 , „ , , 

Note added in this version. In fact, in order to have the agreement, we had to make kde'2 ~ k/2Q instead 
of kde'2 = k/Q as done to compare to other results of the same authors. If the latter factor is used also in this 
case, the coefficient in their Eqn. (3) should be ^ instead of ^ . 
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and 

■ SnkdSMR^ , , , 

<^>= W + (61) 



13.1 Closc-in companions 

In the case of close-in companions, we have £3 — ~e[ — — £g — e'g > 0. We also assume that £5 — 62- 
(The lag £5 corresponds to a radial tidal wave with same frequency as the monthly tide.) Hence 

<e>= -^;^(4 + 144) (62) 

</> = -^^{e'^^^e',) (63) 
For synchronous or stationary companions we have £q < 0(e^) and, thus, 

< ^ > = 2^^'- 

■ ZnSkdMR^ 

2ma^ 2 ^ ^ 

The results obtained are the same as found in the literature since Goldreich (1963), Goldreich and 
Soter (1966) and Peale et al. (1980) if we consider the quality factor Q as the inverse of the phase lag 
of the monthly tide (which coincides with the diurnal tide because of the synchronization J3- 

At the order considered in this paper, the only contribution of the torque due to U22 is a motion 
of the node Ell 

<ON22> = -^^^^^^(\-2e'--S^--e'cos2u). (66) 

However, this precession exists even when 5 — Q. Therefore, it is not a tidal effect due to the misalign- 
ment of U22i but one due to the figure of the body trapped into the 1:1 spin-orbit resonance. 



13.2 Variation in obliquity 

The obliquity / may vary due to the torque acting on the orbit and its counterpart acting on the 
deformed body. Thus, its variation is given by the sum of the two components: < I > and < J >. 
From eqns. pi]) and ([5T|) . we have 



ikdSGA'PR^ 



CQan\ . /_ Cf2an\ , , 
1 - TTT^ )£;,+ ( 1 + ^TT^ ) (eg - £9; 



(67) 



4Ci7a6 

In the case of close-in companions, we have 

CQan C H 

< 1, 

showing that in this case only the part coming from < J > matters. The contribution from < I > can 
be neglected. 

For the conversion of Goldreich and Soter results (as well as for those of Dobbs-Dixon et al. (2004) and 
Mardling and Lin (2004)) we use Q' = 3Q/2fc. 

Since uj is expected to grow monotonically due to the oblateness of the bodies, one part of the given motion 
of the node is not secular. We included it in the given equation for sake of completeness 
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14 Tidal friction in the central body 

Darwin's treatment of tidal friction allowed us to obtain equations that are valid in all cases, provided 
that the tidal phase lags and response factors are left as free parameters. Equations ([56|) . (|60p and 
((6T|) may be used to describe the effects of tidal friction on the central body. However, we have to 
consider separately the case in which the central body is a fast-rotating body (as Jupiter or Saturn) 
or a slow-rotating body (as the Sun and many other main-sequence stars) . 

In the forthcoming equations, one should remember that the parameters R (radius), m (mass), kd 
(dynamical Love number) refer to the deformed body, which is now the central body, while M is the 
mass of the tide-raising body, the close-in companion. The phase lags Si also refer to the central body. 



14.1 Type 1: fl ^ n. Fast-rotating planet 



In this case, the terms appearing in the equations obtained thus far are semi-diurnal (sg, e'j^, £2)1 
diurnal (eg, £g) and monthly (cg) (all positive). The equations for the variation of the elements in this 
case become 



<n> 



3n 
2^ 



< a > 



9n'^kdMR 



5 r 



< e > 



SnekdMR^ 



,27 



(224-34; 



and 



< / > = - 



SnSkdMR^ 



The equations found usually in the literature (e.g. Yoder and Peale, 1981) are 



<n> 



3n 
2a 



< a > 



< e > = 



9n'^kdMR^ 
hlMR^nkd 

g 



(68) 
(69) 

(70) 

(71) 
(72) 



which are obtained from the above equations after neglecting S*^, considering equal the lags and re- 
sponse factors of both diurnal and monthly (radial) tides and introducing |eo| = The coefficients 
in these expressions are slightly different from those found in Kaula(1964) and in Goldreich and Soter 
(1966) since those authors adopt here M -I- to ~ to, which is equivalent to neglect the reaction force. 
There is also a sign difference in Kaula(1964), but it is due to the opposite definition of the corre- 
sponding phases. 



14.2 Type 3: J7 ^ n. Slow-rotating star 

In this case, the terms appearing in the equations obtained thus far are semi-annual (eq, £g), tierce- 
annual {e'l), annual (ej, £5) (all negative except £5) and the "diurnal" lag (£g) (positive). The equations 
for the variation of the elements in this case become 



3n . 9n^kdMR^ r , . 147 , 7 

<n> = -— < a> = - -'^i^^'^' - 

2a 2ma° 



e'{5e'o - —e[ - -e',) (73) 



SnekdMR^ , , 49 , 5 , , 



4TOa5 ^ " 4^4 
where we have introduced £5 = —£3 and £g = e'q. We also have 



■ 3nSkdMR^ 

4TOa5 ^^^^ 
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The equivalent of equations ()7ip and ()72p is obtained assuming that the annual, semi-annual and 
tierce-annual phase lags and Love numbers to be equal. They are 

3n Qn^kdMR^ / 57 

<h> = < d > = ' ' ■ 



2a 2mQa^ 
bMR^n 



<e>= , ^ / e (77) 



where we have introduced |eo| = l/Q- 

One may note that, as expected from the opposite geometry presented by Types 1 and 3, the 
effects on the semi-major axis and eccentricity have different signs. However, they do not have the 
same magnitudes because some phase lags have different signs. The consideration of equal signs for all 
in Type 3 would be physically equivalent to put some waves with a phase advance with respect to 
the corresponding tidal waves in the static model (instead of a lag behind them). 

The above results are only valid if J? 3> n or i7 <C n, respectively. For all intermediary cases, the 
general equations should be used. 

14.3 Rotation of the central body 

From eqns. ([50]) and (pij) . introducing the lags in the same way as above, we obtain: 
14-3.1 Type 1 central body 



<Q> = - 



ikdGhPR^ ( ,^ 15 



2Ca' 



6 



3kdGAPR^ 



tj i\i /] J. r J. J. I, . , , 



14.3.2 Type 3 central body 



<^>= - tele + e^(-20e- + 49e; + 4) + 2S\-e', + 4)) (80) 

■ 'ikdSGhPR^ 

ACna^ (2^0 -£9) (81) 

In the circular approximation, the results given by Goldreich and Soter (1966) for < Q > coincide 
with those given by eqns. ([55)1 and ([75]) if we adopt the same definitions for the moment of inertia used 
by Goldreich and Soter, i.e., C = amR^ and C = O.AMR^, respectively. 



15 Linear theories with a constant time lag 

Several of the existing tidal friction theories introduce ab initio the tidal lag by assuming that the 
tides correspond to the position of the tide-raising body M at a time r before the current time (see 
e.g. Mignard, 1979). For the sake of comparing the results of linear theories to those obtained here, 
we give below the resulting equations when the are assumed to be proportional to the frequencies 
of the corresponding tidal waves, with r as the coefficient of proportionality^. 

For the equations giving the variation of the rotational state of the deformed body, we obtain 



<f2> 



inkdGM^R 



2n5 r 



Ca^ 



{l^'i^^-\s^)-{l+'4^^-\s^)^]r (82) 



Note added in this version. This is equivalent to assume the lags proportional to the frequencies and the 
same dynamical response factors kd for all terms. 
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< J> = 



SnkdGM'^R^ 



5(1 



1 [2, 



-)r. 



(83) 



Cna^ 2n 

Eqn. ([5^ reproduces eqn. (18) of Dobbs-Dixon et al. (2004) to the adopted order in eccentricity, 
except for modifications done there using Kepler's third law, the effects due to stellar wind and some 
specific considerations on the parts of the body participating in the angular momentum exchange . 
The so-called quality factor used there is related to the parameters used here through Q = l/nr. 
For the equations giving the variation of the orbital elements, we obtain 



< n> 



3n 

'2a 



< a> 



5 r 



fl + 23e2 



< e > 



27n'^ekdMR^ 



(1 



<I>= -- 



SnkdSMR^ 



2ma^ 



1 



11 f2. 



[2t. 
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1 



n 

n 



(84) 
(85) 
(86) 



These equations are the same found in several papers on tidal friction on close-in exoplanets using 
Hut's approach (e.g. Mardling and Lin, 2004, Dobbs-Dixon et al. 2004). It is worth noting that the 
results in Hut (1981) are not given by expansions but by closed formulas. 



16 Cumulative orbital variations due to tides in both bodies 

In this section, we add the variations of the mean-motions and eccentricities due to the tides raised in 
both the central body and the companion. In the equations giving the variations due to the tides in 
the central body, we make the substitutions kd = kdA, ^'j — s'jA^ = M = uib, R — Ra, S — Sa 
and J7 = J7^. In the equations giving the variations due to the tides in the synchronous or stationary 
companion, we make the substitutions kd — kds, s'j = s'jb^ m = tub, M = uia, R = Rb, S = Sb and 

n = Hb- 



For sake of simplicity, we introduce the factor 

kdB 



kdA 



2B 



~-0A 



Vms/ \Ra^ 



(87) 



This is the same factor introduced in Yoder and Peale (1981) if we assume that kdA and kdB are 
the Love numbers associated to the tidal waves whose frequencies are 2/2^ — 2n and 2Qb ~ n and 
introduce the ratio of the phase lags instead of the ratio of the quality factors. The conversion is done 
using \e2B/£oA\ = Qa/Qb- 

Let us initially add the variations of the mean-motions and eccentricities given by eqns. (j68p . (|57p . 
([5^ and ([M]) corresponding to a fast-rotating central planet and a satellite in stationary rotation. 

If we assume that the lags are independent of the frequencies and equal, we obtain 



^^^^ _ 9n^kdAmBR\e',^ A ^ 51 _ ^ ^ 

2mAa^ V 4 ^ -"V 



9n^kdAmB R\e'aA , ^ 
' V ^ 4 

■dAnT-BR^SQy^ , 
2mAa^ 4 



2mAa° 

<e>^ ^^-kdAmBR\e',^ ^gg^ 



Inserting kdAs'oA = kA/QA and 5_b = into eqns. ([88]) and ((89)) . they become the ones often 
appearing in the literature on satellites' tidal friction and are the same used by Peale and co-workers 
(e.g. Yoder and Peale, 1981). 

In the case of linear theories with a constant time lag, phase lags are proportional to the corre- 
sponding tidal wave frequencies and the equations for the variation of the elements become 

< e > = 3nefc.^m^^ie-^ 11 _ 

2mAa^ ^ 2 ' ^ ' 
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It is worth noting that there is some robustness in the last two sets of equations due to the fact that 
for fast-rotating central bodies Eq = e[ — e'^. The only difference between the two sets conies from 
the contribution of the monthly (radial) tide (£5) in the central body and the non- vanishing term in 
which do not appear in the first set of equations because of the assumption of the same phase lag 
for the diurnal and semi-diurnal tidal waves. Let us now do the same for a slow rotating star and an 
exoplanet in stationary rotation, but given the diversity of frequencies in this case, we consider only 
the case of tidal lags proportional to the corresponding tide frequencies. We expand eqns. ((84)l and 
([55)) . to the order OiflAjn) (^ 1) and assume Qb — ^{1 + 6e^). Hence (correcting the typos of the 
previous version): 

< n > = ^^^k,AmBRpeU\ / ^ ^ ^^^2 ^ ^(7^2 _^ 3%) + ^ {19e^ + S^)) (92) 

2mAa^ V 9 18 n / ^ ' 

In all cases, if the satellite or exoplanet is in synchronous instead of stationary rotation, we have 
the same equations. However, the giant planets of the Solar System do not show a measurable J22 and 
we should not expect a permanent equatorial ellipticity in hot Jupiters. Therefore, hot Jupiters with 
quality factor not too large will have a fast synchronization but will not reach exact synchronization 
before circularization. They will be rather driven to a stationary rotational state. Anyway, in both 
cases the equations giving tidal variations of the orbital elements are the same. The term in 19De^ 
sometimes found in equations for the tidal variation of the semi-major axis and mean-motion of close-in 
exoplanets does not take into account that a synchronous rotation is only possible if an external torque 
counteracts the tidal torque and may be substituted by TDe^ as given above. The only difference is 
that, in the case of synchronous rotations, we will have to consider also the non-tidal effects due to 
the equatorial asymmetry of the companion. 



17 Motions with respect to an invariable plane 

Throughout the previous sections, we have considered the relative positions of the equatorial and 
orbital planes. None of these planes is fixed. In more general situations (as considered in the next 
section), it may become necessary to refer their positions to an inertial frame. In that case, we may 
adopt the invariable plane (normal to the total angular momentum) as reference plane. If all poles 
are represented on the celestial sphere, the equatorial and orbital poles may be considered as three 
points whose "masses" are the corresponding angular momenta and whose gravity center lies on the 
pole (projected in the origin) of the invariable plane. The distances of the poles of the two equators 
to the pole of the orbit are the mutual inclinations (obliquities) /. The angle between the arcs joining 
the pole of the orbit to the pole of the two equators is the distance of the nodes of the two equators on 
the orbital plane. These data are sufficient to completely determine all parameters of the considered 
planes. 




Fig. 7 Projection of the orbital (Po) and equatorial (Pi,P2) poles on the invariable plane. 5N is the angle 
between the nodal lines. Ii and I2 are the two obliquities. 



We have to consider two steps. The first one is to determine the positions of the orbital (Pq) and 
equatorial (Pi, P2) poles with respect to the pole (O) of the invariable plane when we know 6N, /i, J2. 
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If Mq, Ml, M2 are the angular momenta, respectively, orbital and equatorial, and if we denote by 
Jq^JitJi the distances (angles) of the corresponding poles to the pole of the invariable plane, we 
have: 

Ml sin PiOPo - M2 sin P2OP0 = (94) 
Ml cosPiOPo + Ml COSP2OP0 + ^/o = 



cos/i = COS Jo COS J7i + sin Jo sin J7i cosPiOPo (95) 
cos I2 = cos Jo cos Ji + sin Jo sin J2 cos P2OP0 

and 

cos 1\ cos I2 + sin 1\ sin I2 cos — cos Ji cos J2 + sin Ji sin J2 cos P1OP2 (96) 

which are enough to determine the arcs Jo, Ji, J2 and the angles between them. 

The converse step is to determine Ii,l2, SN when the position of the poles with respect to the pole 
of the invariable plane are known. This can be done easily using eqns. (|95l-[96)). 

The possible polar oblateness of the deformed body was not considered in previous sections. Indeed, 
a torque component directed along the nodal line does not change the given results. However, it will 
make the nodal line precess and such precession (a uniform variation of SN) cannot be neglected when 
considering eqn. ([M)) . 



18 Equations of Motion 

We reproduce below the classical form of the equations of motion of two planets, including the forces 
due to tidal friction. 

Let us consider two bodies with masses tua (primary) and ms (companion). 

Let fo be the attractive force acting on ms due to uia in a point mass model and — fo its reaction 
acting on uia] 

Let fA be the tidal force acting on ttia due to the tidal deformation of m^. It is given by eqns. 
making the substitutions kd = kds, M = m^, R = Rb, S = Sb (sine of the obhquity of the 
equator of ms), s'j = s'jB ^ = (argument of the periapsis reckoned from the intersection of 
the orbit and the equatorial plane of B). Let — be the corresponding reaction acting on m^. 

Let fs be the tidal force acting on uib due to the tidal deformation of ttia- It is given by eqns. (|23l- 
[25]) making the substitutions kd = kdA, M = itib, R — Ra, S = Sa (sine of the obliquity - inclination 
over the orbital plane - of the equator of toa), £j = e'jA^ lo = loa (argument of the periapsis reckoned 
from the intersection of the orbit and the equatorial plane of A). Let — fs be the corresponding reaction 
acting on niA- 

The equations of the motion with respect to an inertial reference frame are: 

ruAVA = -fo + fA-fs (97) 
mBVB = fo - f4 + fs- 

Let, now, r = — be the radius vector of the companion in a reference frame fixed in the 
primary. From the above equations we obtain 

rriBV = (l + !I^)(fo-fA + fi3). (98) 
niA 

The only critical step comes from the fact that the forces f^ and fs were given in different reference 
systems. In previous sections, all results corresponded to scalar quantities - work, semi-major axis, 
eccentricity, obliquities - which are not affected by the frames being used, but, in the above equation, 
it is necessary to take this into account. If [^a] and [f^] represent the vectors whose components are 
given by eqns. ([23]-[25|), we have to change eqn. ((98| into 

ms[r] = (1 + — )([fo] - + M)- (99) 
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where e*'^ indicates a 7r-rotation around the axis normal to the orbital plane. For all other vectors, 
[..] represent their components in the system with center in tua- Note that because of the adopted 
spherical coordinates frames, we have e''^[a, b, c]'^ = [—a, b, —c]'^- 

The other point to be taken into consideration is the relationship between the angle variables in 
both cases. The mean anomaly £ does not depend on the adopted frame, but the arguments of the 
periapsis: tOk are different and their difference is the angle 6N introduced in the previous section. 

18.1 Disturbing Function in perturbation equations 

Many theories since Darwin's (e.g. Kaula, 1964, MacDonald, 1964) use the classical Gauss or Lagrange 
equations for the variation of the orbital elements. In order to assess the consequences here of eqn. 
(1981) . let us consider, for simplicity, the particular case in which the companion of mass ttib is a mass 
point that cannot be deformed. In such case, the equations of motion are simply reduced to uib^ = 
(1 + ^^)(fo + fB). The central force fo derives from the two-body potential Uq and the perturbation 
derives from a potential TZ, i.e. = — ms gradT^. The disturbing potential TZ appearing in Gauss and 
Lagrange equations corresponds to an "external" perturbation acting on ms. However, in the present 
case, the perturbation acting on niB is "internal" to the system of bodies. Therefore, as is usually done 
in the formulation of an N-body problem, the reactions must also be taken into account, that is, TZ in 
Gauss and Lagrange equations has to be substituted by (1 + ^^)TZ to take into account the reaction 
on ruA of its tidal action on m^. This correction can be neglected only when niB ^ m-A, which is the 
case in the cited theories. 

For sake of completeness, we recall that the reaction to fg is duly taken into account in the current 
formulation of Gauss and Lagrange equations. 

19 Conclusions 

We have done a review of Darwin's theory for bodily tides and compared its results to those published in 
many papers dealing with tidal evolution of orbital elements and obliquities and the energy dissipation 
inside bodies due to tidal friction. As expected, our main results coincide with the general results 
given by Kaula (1964) to the order of the approximations adopted. The only differences are that 
we present the results through explicit formulas and that we have not introduced approximations as 
niA + TTiB — TTiA', the given equations remain valid even in the case of larger mass ratios as in the 
Earth-Moon system and in the possible case of a close binary system formed by a normal star and a 
big brown dwarf. 

The comparison to more recent results allowed us to know, in each case, the additional hypotheses 
made by the authors. In some cases the formulas used (e.g. Yoder and Peale, 1981) correspond to 
assuming tidal lags independent of frequencies (as done in eg ns. [551 and [55)1 while others (e.g. Mardling 
and Lin, 2004) correspond to adopting the same linear relationship between lags and frequencies intro- 
duced by Darwin (as done in Sec. [15]). Other consequences are the differences between the equations for 
the variation of the orbital elements in the cases planet-satellite and star-planet, due to the different 
rotational state of the central body in these problems. 

The agreement between the various theories is limited to terms corresponding to the circular ap- 
proximation. Many disagreements are found in terms depending on the eccentricitj{3 showing that 
it is necessary to improve our knowledge of the physics of the processes leading to tidal lags before 
extending the results to higher orders in eccentricity and inclination. In practical applications, these 
divergences are absorbed in the poorly known quality factor Q and results of evolution studies do not 
diverge in the same proportion. 

Emphasis is given in the paper to the case of companions having reached one of the two possible 
final states of the companion synchronization. If tides are the only source of perturbations in the 
system, perfect synchronization cannot be achieved while the orbit is not completely circularized. 
If a remnant eccentricity exists, the final state is rather non-synchronous, the rotation reaching a 
stationary state with a rotation velocity slightly larger than the orbital mean-motion. In order to 
have true synchronous rotation, it is necessary to provide the system with a torque counteracting the 



Other disagreements axe pointed out by Efroimsky (2008). 
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tidal torque. This additional torque can be provided by a permanent solid-like equatorial asymmetry 
of the companion. Tides drive the figure of the body to a small misalignment and the misalignment 
provides one torque acting in opposition to the tidal torque. The results for the tidal variations on the 
system are the same in the stationary and synchronous cases, but, in the synchronous case, we have 
several non-tidal effects coming from the figure of the body which may change the evolution of the 
system. Some equations are found in the literature using the condition f2 = n, which corresponds to 
the truly synchronous case, but without taking into account the counteracting torques allowing for the 
synchronization. Those equations are incomplete. 

Darwin's harmonic decomposition of the tide allowed us to identify the role played in dissipation 
by different tidal waves. In general, dissipation is dominated by the semi-diurnal tide: the harmonic 
whose period is half the synodic rotation period. However when the rotation is captured in the 1:1 
spin-orbit resonance and the synodic rotation period becomes infinite, dissipation is mainly due to 
the monthly tide: the harmonic whose period is the orbital period, which is now equal to the sidereal 
rotation period. In this case the dissipation becomes proportional to the square of the eccentricity. 
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Appendices 

A Equatorial prolateness of Roche ellipsoids 

Roche ellipsoids are figures of equilibrium of a homogeneous fluid with rotation synchronous to their orbital 
(circular) motion about a spherical companion. The analysis of the equilibrium equations shows, however, that 
the angular orbital velocity of the system and the angular rotation velocity of the body appear separately in 
the equations. Therefore the equilibrium equations remain valid even if there is no spin-orbit synchronization. 
The resulting equations in this more general case are (Jeans, 1929, sec. 206): 



2^GpaHc r j^, - - a'n' = (100) 

Jo 



A{a'^ + x) 



^ ^ ,3 f dx GMh ,2^2 ^ 

J A{b2 + x) 7-3 



3 / dx . GMc^ 



(except for terms of the order of the eccentricity) where p is the density, r is the distance between the two 
centers, i7 is the angular rotation of the body and M is the mass of the tide generating body. The axes of the 
ellipsoid axe a, b, c {c <b < a) and = (a^ -|- a;) (6^ -|- x){c^ + x). O is axi auxiliary function, the definition of 
which is not explicitly required in this context. 

The shape of the ellipsoid equator is obtained by eliminating O between the two first equations: 



^ , r {a''-b'^)xdx GM(2a'' + b^) , 2 ,2^^ 

2nGpabc / , , „^ = hx"^ — - + (« " ^ )^ (101) 

Jo " ~ 



A{a^ + x){b'^ + x) 



If the oblatcncss of the ellipsoid is small, the integrals may be evaluated by means of some trivial expansions 
(Tisserand, 1891) to give, to the first order of approximation in the equatorial prolateness: 

^ , n (a^-b'^)xdx AGm 
^'^^^"^'^X Aia^+xW + x) = -m'- 

{R is the mean radius and m is the mass of the ellipsoid). At the same order of approximation, we have 
{a? — b^) = 2E?e. Wc also adopt the approximation 2a^ + fo'^ = (the corrections of the order of the polar 
oblateness may be discarded because of the multiphcation of this term by the small quantity 1/r^). Hence 
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B Variations in the argument of periapsis and nod^ll 



For sake of completeness, we should also consider the variations in the argument of periapsis and the time of 
periapsis. They can be easily obtained from the corresponding Gauss equations (see Beutler, 2001, p. 230). The 
force components Fr, Fs, Fw (radial, transverse and normal)(3 appearing in Gauss equations are related to 
the forces Fi, F2, F3 through the same rotation shown in fig. I5h°l 



Fr = Fi 

Fs = -F2 sin 13 + F3 cos /? 
Fw = —F2 cos P — F3 sin /? 



(104) 



We introduce these relations in the corresponding Gauss equation multiplied by (1 + M/m) as discussed in sec- 
tion ll8.l1 and average the result over an orbital period. We note that, at variance with the other perturbations, 
in the case of the argument of perihelion we have to consider both the dynamic and the static tide. Hence 



QkduMR^ 



(3eo - £2 " £5 + 3£6 " Se's - e'15 + e'le + 3e25) 
- -^(Geo - £2 - £5 - 6£7 - 12£8 + £'11 + £12 - 2e'i5 + '^e'la) 



sin 2uj 



and 



2maB o^,. (105) 

respectively. The comparison of them shows that the part coming from the dynamic tide is attenuated by the 
product by the lags e'i and by and 5*^ so that the motion of the periapsis is dominated by the static tide. 
The new lags appearing in the above expressions are given below: 



Subscript No. 


Frequency 


11 


20 -n 


12 


n 


14 


n + n 


15, 16 


Q — n 


17 


Q — 'in 
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n -2n 



The quotient of the contribution from the tide raised in a close-in companion to that coming from the tide 
raised in the central body is 

kdB fPA\^RA 

kdA ^Pb Rb 

where A, B represent the central body and the companion respectively and p are the mean densities. It shows 
that the effects of the tides raised on the companion dominate < uo >. 

The comparison of the tidal motion of the periapsis and the mean- motion is roughly given by < uj^jn >. 
It is of the order 0{R%/a^). (For this estimate we assumed similar densities and Ra ~ a.) 

In the case of one companion in synchronous rotation, we have to consider also the contribution coming 
from the asymmetry of the companion. If we proceed in the same way as above, but using U22 instead of U2, 
we obtain 

3nR% J22 



< UJ22 > ~ 



4a2 



(I6 - 15e^ - + 9(e^ - 5*^) cos 2lu) 



(106) 



whose main part is again of the order 0(R%/a^). However, this result is given only to show the order of 
magnitude of the results. The main contribution from the shape of the companion is expected to come from 
the zonal term proportional to J2P2{cosd)/r^ , which has not been included in the calculations presented in 
this report. 

With the corrections published in the Errata. The results given in the paper are incomplete. To obtain 
complete second-degree results in e, /, it is necessary to consider in U2 terms up to fourth degree. In particular, 
it is necessary to use 



V + u> — ^ sin(2t; + 2uj) sin^ / — sin(2(i; + 2d) — ^ sin(4tt) + 4v)j sin* /. 



instead of the approximation given in eqn. (8), as well as higher-order equations for the Keplerian motion. 

The higher-order terms thus introduced also contribute the variations in the node given in the added sub- 
section. 

Usually these equations are written in terms of the accelerations R, S, W instead of the forces 
The second- and third-degree contributions to Fi,F2,F3 are given in section [6] and Appendix D of the 
astro-ph version of the paper. 
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B.l Variations in the Node 

At the same approximation, the motion of the node is given by 
— 9kdnMR^ 



< ONlag > = 



16ma^ 



(3eo - £2 - £5 + See - Ses - £15 + e'le + ^£25) e sin 2uj. 



(107) 



Because of the rotational symmetry around the axis a (see Fig. [TJ, the static tide does not contribute a 
variation of the node. 



C Forces arising from (7iag0 



9kdGMM*R^ 

8^3^*4 



P^e'o (2 - 5e^ - S^'j sm(2if* - 2£ - 2uj) 



+eP^(l - ^5"^) {le'^ sin(2<^* - 3f - 2lj) - e'2 sin(2vp* - I- 2uj)j 

+ 17e^P^e3 sin(2<^* - 4^ - 2u) + P'^S'^e'i sin 2ip* 
-ee5(4 - 6P^ + 3P^S^) sin£ - 3e^£6(2 - 3P^) sin2£ 

+P^S'^e'7 sin(2£ + 2uj) + 2Qs{e'^ COs(<^* — 21 — 2lo) — £g cos if 
+ ieP^S'^ (^3e'io sin(2v3* + I) + 3eii sin(2<^* - tj 

-e'la sin(^ + 2lo) + Te'va sin(3£ + 2lo)^ 
—Qesi^e'n cos{(p* + £) + 3e'i5 cos{(p* ~ £) 

+£'ie cos{(p* — £ — 2uj) — 7e'i7 cos(v5* ~ 3£ — 2lu)^ 



(108) 



3kaGMM*R^ 



ga3j..4 



'Eo (2 - 5e^^ sin(2<^* - 2^ - 2cj) 

+ ieQ(^7ei sin(2(^* - 3£ - 2lj) - £3 sin(2<^* - £ - 2uj) 
+ 17e^Q£3 sin(2v3* - 4f - 2uj) + GeQe'r, sml + de^Qe'e sin 2£ 
+4(1 - 2P^)S(^£8 cos(<y3* - 2£-2uj) - £9 cosi^*^ 

-2(1 - 2P^)es(se'i4, cos((^* + £) + Se'15 cos{(p* - £) 

+£i(i cos((/3* — £ — 2lj) — 7£'i7 cos{(f* — 3£ — 2lo) 



(109) 



?,kiGMM*R^ 



£0 (^2 - 5e^ - 5"^^ cos(2v3* - 2£ - 2uj) 

+2eP(l - is^) (^7£i cos(2v7* - 3^ - 2lj) - e'2 cos(2(^* - £ - 2a;) j 
+34eV£3 cos(2<^* - 4£ - 2tj) + 2PS^£l cos 2^?* 



+eP5'' 

+ %eSi 3£'i4 sin(</?* + ^) + Se'is sin(v3* — £) 



^£8 sin((j9* ~ 2£ — 2lu) — £9 sin tp*^ 
' (^3£io cos(2v9* + ^) + 3£ii cos(2v9* - £)^ 



+£i6 sm{(p* — £ — 2lj) — 7e'n sm{(p* — 3£ — 2Lo)j 



(110) 



The appendices C, D, E were published only in the astro-ph version of the paper 
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We have used here ^ = 2 — AP^ and = Q. The latest derivative leads to the formation of several terms 

proportional to QS^, which were dropped out from the expression of Fj''^ (but included in the third-degree 
expressions of Appendix D). 



D Third-degree forces 

The terms of third and fourth degree of f/2 are: 



S3U2 



SkfGMR^ 



ie^f 18cos^ + cos3£ + cos(2vJ* - £ - 2oj) - 123cos(2(p* - 3£ - 2a;) 
8 V o 



+ ^ cos(2(/3* +£-2ui) + ^ cos{2(fi* - 5£ - 2ui)) 
+^eS'^(^- 7cos(2<^* - 3£ - 2w) + cos(2(p* - £ - 2a;) - 6 cos £ 

+ 3 cos(293* +£) + 3 cos(2v?* - £) - cos(£ + 2a;) + 7 cos(3£ + 2a;)^ 

+eQs(^sm{ip* - £ - 2a;) + 3sin(<^* + £) +3sin(v5* - £) - 7sin(</j* - 3£ - 2a))^ 



(111) 



and 
54U2 



8a3r*3 



^e''(l5-|-28cos2£ + 77cos4£) + ie''(l3cos(2(p* - 2£ - 2a;) (112) 
- — cos(2(p* - 4£ - 2a;) + 533 cos(2vj* - 6£ - 2a;) + - cos(2(p* + 2£ - 2a;)) 

+ ie^5'^ ( - 6 - 18 cos 2£ - 10 cos(2£ + 2a;) + 34 cos(4£ + 2a;) + 6 cos 2ip* 

+ 10cos(2<^* - 2£ - 2a;) + 9cos(2(y^* + 2£) + 9cos(2¥;* - 2£) - 34cos(2¥p* - 4£ - 2a;)) 
+ iQSe^(6sinv?* + 9sin(v5* + 2t) + 9sin((p* - 2€) + 10sin(v5* - 2£ - 2a;) - 34sin(¥;* - 4£ - 2a;)) 
+ iQS'^ ( - 2 sin </p* + sm{ip' + 21 + 2uj) + sin^Lp' -21- 2a;)) 
+is"' ( cos(2(^* + 2£ + 2a)) - cos(2(^* - 2£ - 2a))) + ^Q^S^ (l - cos(2£ + 2a;)) 

where we have adopted the approximation = sin^ 0* ~ 1 when this quantity appears multiplied by a third 
or fourth-degree term, because (1 — P^) = cos^ Q* is of the order of . 
The terms of third and fourth degree of C/iag are: 



8a3r*3 



^ e^(27e5 sin£ -|- 53e'i8 sin3£) + sin(2v?* - £ - 2a;) 
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(113) 



123£'i sin(2¥;* - 3£ - 2a;) -|- ^eig sin(2¥;* - 5£ - 2a;) -|- ^£20 sin(2(^* +1- 2a;)) 



+ -e5^ ( - 665 siu£ - e'la sin(£ + 2<J) + 7ei3 sin(3£ + 2uj) + 3eio sin(293* + i) + Se'n sin(2vp' - i) 

+ e'2 sin(299* - ^ - 2a)) - 7£i sin(2v3* - U~ 2a))) 
— SQe^e'ig cos((p* — £ — 2a)) — 7e'i7 003(95* — 3£ — 2a;) + 3e'i4 cos((p* + £) + Se'is cos(v9* — £)) 



and 



3kdGMR^ 



8a^r 



7 11 1 "^Q 

^e\e'e sin2£ + ^^£^2 sin4£) + ±e''{^s'o sm{2<p* - 2£ - 2a;) 

1599 , 



(114) 



-I- £31 sin(2(^* -I- 2£ - 2a;) - 460£3 sin(2¥;* - 4£ - 2a;) + e'29 sin(2¥;* -61- 2a;)) 

- ISeg sin 2£ - lOey sin(2£ -I- 2a)) + 664 sin 2ifi* + 10e'„ sin(2(p* - 2£ - 2a;) 
+ 34£2i sin(4£ + 2ui) + 9e27 sin(2<^* + 2£) + 9e28 sin(2<^* - 2£) - 34e3 sin(2(p* - 4£ - 2a;)) 
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+ ^QSe^ ( ^ 6eg cos ip* — 9e24 cos{ip* + 21) — 9e25 cos(</7* — 2(,) + 34e23 cos((y3* — 4^ — 2tt;) 

- lOeg cos(v9* - 2^ - 2i^)) + igS^ (269 cos p* - e'ae cos(<^* + 21 + 2u;) - eg cos((^* -21- 2uj)) 

+ is*(e3o sm(2</p* + 2£ + 2i^) - ed sin(2v9* - 21 - 2iu)) - ie7Q^S%m(2^ + 2tj) 
8 2 

In these expansion several new lags Sj appear. They are not individually identified (except for those given 
in Appendix B). By definiton, they are the lags of the arguments of the terms in which they appear (see Sect. 
4). For instance, in the term [...]e'ii sin(2(/p* — i), en is the lag of the argument 2^5* — £. 

The third-degree components of the forces are 



^ e ( - Sleo + —El + -^-£2 - 68e3 + -es + Gee j sm I 

+ eS^ ( — 8ed + 7ei + £2 — 6e5 + 8eg + 3ei4 — 3£i5 — e'la — Te'i 

^ 3/143 , 7 , 43 , ^„ , ^„ , ^ 53 , 845 , ^ 1 
+ e (—£0 - j£i - —£2 + 7£5 + 6e6 + — £13 - —£19 + - 

+ eA I - 2eo + -£2 - -£5 + 2£7 + 4eg + -en - -£12 - 3ei5 - 
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'17) sin£ 

43 , „ , „ , 53 , 845 ,. 1 

3 4""^ ^£2 + 7£5 + b£6 + — ^18 24~^^ 



24'^'" ^ 2A^''^°^ ^^^^^ 



' ^ 3 " 4^ 4 " ^ ° 12 '° 24 24 
+ eS [ - 2eo + -£2 - -£5 + 2£7 + 4eg + -en - -£12 - 3ei5 - £1, 

eS^ ( - 2ed - ^e'l + 4e4 + ^eg + 2e7 + 4e8 - 8e9 + |eio + ^eis - 3ei4 



■'le) sin(^ + 2lj) 

' 7e'i7) sin(3£ + 2uj) 



S3F2 = 



26^5* (Sed — 7£'i — £2 + 6e5 — 4eg + 4e9 + 3ei4 + 9e'i5 + e'lg — 21e'i7) cos(^ + to) 
+ e'^S(l3e'o + Te'i - 3£2 — 17£3 + 6e'r, + Qe'e — ^eg + -^eg + Geig + Geie) cos(^ — ui) 



(116) 



1 1 CO 

e^5'( + 29e;) - 7e'i - 5£2 - 17£3 + 18£5 + 9£6 + Y£g - y£9 - 18£i4 + 14e'i7) cos(3^ + uj) 
(e[) — e4 — e7 — |eg + ^£9) ( cos(^ + lj) — cos(3/ + 3lu)) 



+ 



8a7 



9e^5'£25 cos(^ ~Lu) + 2e^S(3£9 + 5£8) cos(£ + uj) + S^(2e'-j + £» - 2eg) cos(^ + to) 



^S'(34£'23 - 9£24) cos(3^ + Lo) - {2e'j - £26) cos(3£ + 3a;) 



.53-^3 



SkdGM'R-' 



_ , 423 , 7 , 
( + 58£o — £1 - -£2 



136£3) cos^ 



(117) 
(118) 



+ eS^ {^2e'o + ^£1 ^ ^£2 — 8£g — 8eg — 3ei4 — Se'i^ + e'm — 7£i7) cos^ 

3/_La8 ' _L ' 43 , ^ 845 , ^ 1 , ^ „. 
- e ( + 98£o + -£i - y£2 + -[^^19 + Y2^2o j cos 3£ 

+ eS^ ( — 4ed + ^e'l + ^e2 + 4eg — 3£ii + 3£'i5 — £'16) cos(^ + 2uj) 
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+ eS ( — —e'l — -£2 — 8£4 + 4e8 + 8£g — 3£'io + 3£'i4 + 7e'i7) cos(3^ + 2lu) 

In the expression of 53F2, one may note two parts graphically separated. The first of them comes from 
SsUiag and the seco nd from (54L'iag. In the same way, 53f7iag contributes terms of degree 2 to F2, which were 
included in eqn. 1241 

The lost of one degree in the derivation of the force comes from derivatives with respect to Q (Q is of the 
order 0{S), but its derivative 2 — 4P^ ~ —2 is a finite quantity). 



E Coriolis accelerations 

In Darwin's theory, the rotation of the tidally deformed body does not appear explicitly in the tidal forces. 
It only enters through the tidal lag angles and by the association of each lag to the corresponding wave in 
the body. For instance, when we say that e2 corresponds to the semi-diurnal wave (in Type 1 tides), we axe 
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implicitly assuming the equator of the deformed body as reference plane. Since this plane is not fixed, we have 
to take into account the rising of apparent forces due to its motion (see e.g. Ferraz-Mello, 2007, p. 17): 

Fapp = AfO X (r X O) + 2Afv x O (119) 

where O is the rotation angular velocity vector of the reference frame The first term is the centrifugal force. It 
is proportional to the square of the angular rotation velocity and can be neglected at the order considered in 
this report. Indeed, in this case, the frame rotates around the nodal line with |0| = J ~ 0{Se). Therefore the 
centrifugal force is of order 0{S'^e^). The second term is the Coriolis force. This term is of the same order of 
the terms considered in previous sections and, for sake of completeness, needs to be calculated. 

The vectors v and O belong to the orbital plane (O lies on the nodal line). Therefore the Coriolis force 
is perpendicular to the orbital plane. The first conclusion is that no work is done (the scalar product by v 
is zero) and, thus, it does not affect the orbital semi-major axis. However, it has a torque directed along the 
transverse component of the orbital velocity. First we have to express v in a reference system whose in-plane 
unit vectors are n along the nodal line and ih normal to it. The frame angular velocity vector is then O = Jfi 
and the orbital velocity of M is 



: sin((x; + v) + e sin a;) n — ( cos(lj + v) + e cos a;) m] . (120) 



Hence 

fcor = 2Mv X O = 2Mj^^^=(cos(tj -I- v) + ecosu) u (121) 
Vl — ^ 

where, as before, u is a unit vector normal to the orbital plane. The torque due to the Coriolis force is 

Cor = r X fcor = 2Mj—^^=(cos{Lj + v) + ecoscj) t (122) 
\/l — e2 ^ 

where t is the transvers^3 unit vector: 

t — e**-"^"-* rh = — sm{LO + v)n + cos{lo + v) rh. 
If the two components of the torque are averaged over one period, we obtain, respectively: 

M„ = i MJna^e^ sin 2tj (123) 

Mm = M.jna^{l -I- ^e^) + ^Mjna^e'^ cos 2uj 

To remain in the order of magnitude of the previous calculations, we neglect terms in Je^. Thus M Cor = 
MJna? m. If this torque is substituted in the equation = there results a homogeneous equation in 
the variables J^l-C], N, {MJna^ + i\C\), where A'^ — ON, with the trivial solution 0,0,0. Then |£| remains 

constant, that is, the eccentricity is not affected by the Coriolis force (at the considered order), = (no 
precession of the node) and 

■ M.Jna^ ■ 

ICor = J (124) 

This is the same result which would be obtained from an elementary geometric reasoning. The decrease of the 
obliquity / is equal to the amount of the rotation of the reference plane. In other words, the orbital plane is 
not affected. 

The effect of the Coriolis force on the rotating mass m is similar. Elementary calculations allow us to obtain 

for the apparent torque, M Cor = —JCfl x n which corresponds to counter-rotate the body around the nodal 
line with the same rotation angular velocity as the frame. 
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